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Abstract. The Pierre Auger Observatory detects the highest energy cosmic rays. Calorimetric measure- 
ments of extensive air showers induced by cosmic rays are performed with a fluorescence detector. Thus, 
one of the main challenges is the atmospheric monitoring, especially for aerosols in suspension in the at- 
__ mosphere. Several methods are described which have been developed to measure the aerosol optical depth 

profile and aerosol phase function, using lasers and other light sources as recorded by the fluorescence detec- 
tor. The origin of atmospheric aerosols traveling through the Auger site is also presented, highlighting the 
effect of surrounding areas to atmospheric properties. In the aim to extend the Pierre Auger Observatory 
to an atmospheric research platform, a discussion about a collaborative project is presented. 

Oh 
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^ 1 Introduction 

in 

The Pierre Auger Observatory [TJ[2] is the largest operating cosmic ray observatory ever built. It was designed to 
OA measure the flux, arrival directions and mass composition of cosmic rays from 10 18 eV (electron-volt) to the very 
highest energies. When cosmic rays enter the atmosphere, they induce extensive air showers composed of secondary 
QO particles. Charged particles excite atmospheric nitrogen molecules, and these molecules then emit fluorescence light 
in the 300 — 400 nm range [2|4]. At the Pierre Auger Observatory, the atmosphere is used as a giant calorimeter, 
O^l representing a detector volume larger than 30 000 km 3 . To minimise as much as possible the systematic errors of the 
fluorescence measurements, atmosphere properties have to be continuously monitored [5] . During the development of an 
extensive air shower, the production rate of fluorescence photons depends on the temperature, pressure and humidity 
of the air [HBIB]. Then, from their production point to the telescope, these photons can be scattered by molecules 
(by Rayleigh scattering) and/or atmospheric aerosols (by Mie scattering). Thus, to track atmospheric parameters, an 
extensive atmospheric monitoring system has been developed that covers the whole array. 

The different experimental facilities and their locations are shown in Fig.[T] Atmospheric properties at ground level 
are provided by a network of five weather stations located at each fluorescence detector (FD) site and at the Central 
Laser Facility (CLF). They furnish atmospheric state variable measurements every five minutes. Also, meteorological 
radio-sonde flights with balloons have been operated (for more details, see [10]). For the aerosol component, the 
central laser facility fires 50 vertical shots every 15 minutes during FD operations. Fluorescence telescopes, recording 
the Ultra- Violet (UV) laser tracks, are able to deduce the aerosol optical depth at different altitudes. In November 2008, 
a second laser facility called the eXtreme Laser Facility, or XLF, was deployed. These aerosol profile measurements 
are complemented by cloud measurements by four elastic backscattering lidars located at each eye (for more details, 
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Fig. 1. Atmospheric monitoring map of the Pierre Auger Observatory (from [9j). Gray dots show the positions of 
surface detector stations (SD). Black segments indicate the fields of view of the fluorescence detectors (FD) which are located 
in four sites, called Los Leones (LL), Los Morados (LM), Loma Amarilla (LA) and Coihueco (CO), on the perimeter of the 
surface array. Each FD site hosts several atmospheric monitoring facilities. 



see [H]). A Raman lidar currently under test in Colorado (USA) is scheduled to be moved to the Auger Observatory 
for the Super- Test-Beam project [T^]. To improve our knowledge of photon scattering on aerosols, two Aerosol Phase 
Function monitors (APF) have been installed at the Coihueco and Los Morados FD sites. The APF instruments 
generate a collimated horizontal light beam produced by a Xenon flasher. The light passes in front of one FD site. The 
aerosol attenuation depends on the incident wavelength. This measurement is the main goal of two optical telescopes in 
Auger, the Horizontal Attenuation Monitor, or HAM [13], and the (F/ph)otometric Robotic Telescope for Atmospheric 
Monitoring, or FRAM [H] . 



2 Aerosol effects on light propagation in the atmosphere 

Although the atmosphere is mainly composed of molecules, small particles such as dust or droplets may be in sus- 
pension. These particles are called atmospheric aerosols and their typical size varies from a few nanometres to a few 
micrometres. Most of the atmospheric aerosols are present only in the first few kilometres above the ground. Unlike 
the molecular component, the aerosol population is highly variable in time and location, depending on the wind and 
weather conditions. Aerosol conditions affect the propagation of UV light from air showers. Even if the absolute flu- 
orescence yield remains the largest source of uncertainty for FD measurements [15] . aerosol effects also contribute 
significantly to systematic uncertainties. Absorption in the air is usually negligible for photons having a wavelength 
in the UV domain. Thus, attenuation will represent only the scattering phenomena. The main atmospheric attenua- 
tion processes are Rayleigh scattering, by the molecular component, and Mie scattering, by the aerosol component, 
both being elastic. In the following, only scattering on aerosols will be discussed. For more details on the molecular 
component of the atmosphere, see [10] , 

Two main physical quantities have to be estimated to correct the effect of the aerosols on the number of photons 
detected by the telescopes. These are the aerosol attenuation length as a function of height, linked to the aerosol 
optical depth, and the aerosol scattering phase function. 
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2.1 Aerosol transmission - Mie regime 

Aerosol attenuation of light from an air shower to the FD telescope can be expressed as a transmission coefficient 
r a , or optical transmittance, giving the fraction of incident light at a specified wavelength A that passes through an 
atmosphere of thickness x. If r a is the aerosol optical depth, then r a is estimated using the Beer-Lambert law 

r^ix, A) = exp [—r a (x, A)] , (1) 

under the assumption of horizontal uniformity. Usually, the aerosol population can be described as a superposition of 
several horizontally uniform layers, so the spatial dependence of the aerosol optical depth can be reduced to an altitude 
dependence. Hence, the aerosol transmission from an altitude h to the ground through a slanted path of zenith angle 
9 is 

a(M)" 



r a (M,A) = (1 + H.O.) exp 



cos ( 



(2) 



where H.O. represents a higher-order correction that accounts for the single and multiple scattering of photons into 
the field of view of the telescope. 

Therefore, knowledge of the aerosol transmission parameters at the Pierre Auger Observatory will be acquired 
through measurements of the aerosol optical depth r a (/i, A) at different location in the Auger array, throughout the 
night. Usually, the aerosol concentration decreases rapidly with the altitude. The aerosol optical depths are measured 
in the field at a fixed wavelength Ao- To evaluate the aerosol extinction for a given incident wavelength, a common 
parameterization used is a power law due to Angstrom, 

r a (h,X)=r a (h,X ) 0^Y, (3) 

where 7 is known as the Angstrom coefficient. This exponent depends on the size distribution of the aerosols. When the 
aerosol particle size approaches the size of air molecules, 7 should tend to 4 (mainly dominated by accumulation-mode 
aerosols), and for very large particles, typically larger than 1 /Ltm, it should approach zero (dominated by coarse- 
mode aerosols). Usually, a 7 ~ is characteristic of a desert environment and the aerosol optical depth is more or 
less independent of the wavelength. Measurements at the Pierre Auger Observatory gave 7 values close to zero (see 



Sect. 3.31 



2.2 Angular dependence of aerosol scattering 



Due to the collecting area of the FD telescopes and the typical distance between a telescope and an air shower, only a 
small fraction of the photons produced in the air shower is detected without scattering. Thus, the scattering properties 
of the atmosphere need to be well estimated. The angular dependence of the scattering may be described by a phase 
function P((), defined as the probability per unit solid angle for scattering out of the beam path through an angle 
£. Following the convention for the atmospheric domain, er a -1 dcr a /di7 is the normalised differential aerosol scattering 
cross section, which is identical to the aerosol phase function -P a (£)- The integral of P a (£) over all solid angles has to 
be equal to unity. 

Whereas the molecular component is described analytically by the Rayleigh scattering theory [IB], the Mie scat- 
tering cannot be described by a basic equation for the aerosol component. This is due to the fact that the scattering 
cross section depends on the size distribution and shape of the scattering centres. Typically, forward scattering dom- 
inates in the Mie regime and the forward-backward ratio can vary strongly with aerosol type. More details on the 
relationship between scattering phase function and aerosol size can be found in |17j where a phenomenological ap- 
proach is developed. The aerosol phase function used by the Pierre Auger Collaboration is parameterised by a modified 
Henyey-Greenstein function [TH] 
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4?r 



3 cos 2 C - 1 



(l + . 9 2 -2.gcosC) 3/2 2(l+fJ 2 f 2 



(4) 



where g = ^hg — (cos £} is the asymmetry parameter and / the strength of the second component to the backward 
scattering peak (/ = meaning no additional component). The asymmetry parameter provides the scattered light 
intensity in the forward direction: a larger g means more forward-scattered light (see Fig. [2]) . Values go from g = 1 
(pure forward scattering) to g = — 1 (pure backward scattering), with g — meaning isotropic scattering. The 
second term in the expression - a second-order Legendre polynomial, chosen not to affect the normalization of the 
phase function - is introduced to describe the extra backscattering component. At the Pierre Auger Observatory, the 
goal in monitoring the aerosol phase func tion is to estimate the {<?, /} parameters, two observable quantities depending 
on local aerosol properties (see Sect. 3.2 1. 
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Fig. 2. Henyey-Greenstein functions representing the aerosol phase function for different asymmetry parameters 
gHG and backward factors f. The Rayleigh phase function, proportional to (1 + cos 2 Q an d representing scattering properties 
for the molecular component of the atmosphere, is also plotted. 

3 Aerosol measurements at the Pierre Auger Observatory 

Several instruments are deployed at the Pierre Auger Observatory to observe aerosol scattering properties. The aerosol 
optical depth is estimated using UV laser measurements from the central lasers and scanning lidars (more details on 
lidar measurements in the aerosol phase function is determined with Aerosol Phase Function monitors and the 

wavelength dependence of the aerosol optical depth is measured with the optical telescopes HAM and FRAM. In 
addition, an aerosol sampling program based on filters was undertaken from June to November 2008. 



3.1 Aerosol optical depth measurements 

At the Pierre Auger Observatory, the aerosol optical depth is measured through the night using the fluorescence 
detector to measure light from the central laser facility, located on-site towards the centre of the SD. The main role 
of the CLF is to produce calibrated laser "test beams" . It is powered by a battery, charged by solar panels [Tijj . The 
main component is a laser with a wavelength fixed at 355 nm, in the middle of the nitrogen fluorescence spectrum 
produced by air showers [3]. The pulse width of the beam is 7 ns and a maximum energy per pulse is around 7 mJ. 
This is of the order of fluorescence light produced by a shower with an energy of 10 20 eV. To estimate the relative 
energy of each laser pulse, a portion of the beam is diverted into a photo-diode detector. The analysis of the vertical 
aerosol optical depth uses a vertical beam. A steering head mounted on the roof of the container makes it possible to 
direct the beam towards any direction above the horizon. 

When a laser shot is fired, the FD telescopes collect a small fraction of the light scattered out of the laser beam. 
The recorded signal is not constant but depends on the atmosphere properties. Thus, a method has been developed 
in the Auger Collaboration to estimate the vertical aerosol optical depth T a (h, Ao) with the CLF, where Ao is the CLF 
wavelength [20]. The method to determine r a (h, Ao) normalises the measurement of the laser beam to the signal that 
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C = f + a 




L = 25 - 40 km 

Fig. 3. Geometrical arrangement, viewed from the side, of the central laser CLF and the FD telescope. The 

light is scattered out of the laser beam at a height h corresponding to an elevation angle a and a scattering angle £ = tv/2 + a. 
(_T m _T a )i and (_T m _T a )2 are the total attenuations from the CLF to the scattering location and from the scattering location to 
the FD, respectively. The transition between the Mie domination and the Rayleigh domination is located just a few kilometres 
Above Ground Level (AGL), dependent on the aerosol conditions. 



would be recorded under aerosol- free atmospheric conditions. The molecular references are average CLF laser profiles 
measured during very clear nights, the so-called "reference clear nights" . Once a light profile, averaged over one hour, 
is normalised using a clear-night reference, the attenuation of the remaining light is due to aerosol scattering along 
the path from the laser beam to the FD telescope. A horizontal uniformity for the molecular and aerosol components 
is assumed. The amount of light from the laser beam reaching the detector at the elevation angle a is written as 

(5) 



N obs {a) = N (P m r a )j [P m (I + + Pa (f + «)] ( r m Pa) 2 , 

where N is the number of photons produced per laser pulse, {P m ,P a } are the molecular and aerosol transmission 
factors, and {P m ,P a } are the molecular and aerosol phase functions (see Fig. [3]). The indices 1 and 2 correspond to 
the way from the CLF to the scattering location, and from the scattering location to the detector. 

In the case of a clear night, the aerosol transmission factor P a is equal to one and the scattering over aerosols is 
negligible. Eq. ([5| is reduced to molecular part only 



P, 



(r m ) 2 



(6) 



Applying the aerosol horizontal uniformity condition, P a (/i, a, Ao) = exp [— r a (/i, Ao)/ sin a], and combining A?obs 
and A^noi, Eq. ([5| becomes 



i{h, A ) = - 



1 + sin a 



In 



^Vobs(q) 

Nmol(a) 



-hi 



P a (f + o) 
P m (f + a) 



(7) 



In the scattering angle range seen by the FD telescopes (between 90° and 120° here) the aerosol scattering contri- 
bution is much lower than the molecular contribution. Eq. (Fn) simplifies to 



sin a fN nml (a) 
T&{h, A ) = : In 1 



1 + sin a 



Nobs (a) 



(8) 
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Fig. 4. Aerosol optical depth at 3.5 km above the fluorescence telescopes at Los Leones, Los Morados and 
Coihueco. The wavelength of the CLF laser is 355 nm. Data were acquired between January 2004 and December 2010. 



With these approximations, the aerosol optical depth formula depends only on the elevation angle a of each laser 
track segment, linked to the altitude by h = L tana + /lAugcn where L is the horizontal distance between the FD site 
and the CLF, and ZiAuger the altitude of the Auger array above sea level. 

Figure [4] depicts the aerosol optical depth distribution recorded at Los Leones, Los Morados and Coihueco between 
2004 and 2010. Measurements from Loma Amarilla are not available due to its large distance from the CLF site. The 
XLF, closer to Loma Amarilla, provides optical depth information for this site. Typically, at 3.5 km above the FD 
level, the mean value for r a (ft., Ao) is around 0.04. Note that a cut at 0.1 is used in the air shower reconstruction as a 
quality cut: all the events occurring in a night with a T a (h, Ao) > 0.1 are rejected. 

Horizontal uniformity is usually assumed for the molecular component. Figure [5] shows the scatter plots of the 
aerosol optical depth T a (/i, Ao) measured at Los Morados and Coihueco with respect to the ones measured at Los 
Leones. The altitude is fixed at 3 km above the FD level and the wavelength at 355 nm (CLF wavelength). This means 
that different altitudes are probed at Los Leones and Coihueco, the latter being almost 300 m higher. The optical 
depth data for the different sites suggest that aerosol conditions differ with location. The agreement is better between 
Los Leones and Los Morados, two sites differing only slightly in altitude. The difference could be also explained by 
a different composition of the soil between the Coihueco site and the other ones, aerosols coming mainly from the 
ground. To take into account this non-uniformity, the implementation of the aerosol parameters into the Auger Offline 
software divides the Auger array into five zones centred on the midpoints between the FD buildings and the CLF. 
Horizontal uniformity of the aerosol component is assumed in each slice. Then, each region is divided vertically into 
layers, each layer having a thickness equal to 200 metres. 



3.2 Aerosol scattering measurements 



The FD reconstruction of the cosmic ray energy has to take into account not only the light attenuated during propaga- 
tion, but also has to remove the multiple scattering component adding to the fluorescence light contamination. Aerosol 
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Fig. 5. Horizontal uniformity of the aerosol optical depths measured with CLF shots at Los Leones (LL), Los 
Morados (LM) and Coihueco (CO). Scatter plots show the correlation between two FD sites, the dashed line representing 
full agreement between two data points and the solid line the fit to data (from [5]). 



scattering is described by the aerosol phase function P a (C) which can be parameterised by a modified Henyey-Greenstein 
function as in Eq. @. At the Pierre Auger Observatory, two Aerosol Phase Function monitors, in conjunction with 
the FD telescopes, are used to measure the parameters {<?, /} on an hourly basis during FD data acquisition. The 
APF light sources emit a near-horizontal pulsed light beam in the field of view of their nearby FD site at Coihueco 
and Los Morados. Each APF building contains collimated Xenon flash lamp sources, firing an hourly sequence of 
350 nm and 390 nm shots. The aerosol phase function is then reconstructed from the intensity of the light observed by 
the FD cameras as a function of scattering angle C 7 for angles between 30° and 150°. After corrections for geometry, 
attenuation and collection efficiency for each pixel, the binned APF signal S(() observed is subjected to a 4-parameter 
fit 



s(0 = c 



^m(^Auger) ^at^ Auger) 



(9) 



where {C/ A m (h Auger), C/yl a (/iAuger), 9, /} ar e the fit parameters and -P m (C) the Rayleigh phase function [21.. The 
first two fit parameters can be used to estimate the molecular attenuation length and the aerosol attenuation length, 
respectively; while g and / are used to estimate the aerosol size distribution. 

Figure [6] (left) shows the distribution for the asymmetry parameter measured at Coihueco between June 2006 and 
July 2008. The average g value is (g) = 0.57 ± 0.10, when one excludes the g = nights corresponding to aerosol-free 
nights. Such a value for the asymmetry parameter corresponds to values usually measured in desert locations with 
significant levels of sand and soil dust [22]. The second parameter of the modified Henyey-Greenstein function, the 
backscatter coefficient, is estimated as (/) = 0.40 ±0.10. It means that a backward peak for the aerosol phase function 
exists. Hence, at the Pierre Auger Observatory, the typical aerosol phase function is described by P a (C\g = 0.6, / = 0.4). 

Recently, a new method based on very inclined shots fired by the CLF was developed (laser shots with zenith 
angles higher than 86°) [23]. Following the same idea as previously, knowing the geometry of the laser shot and the 
signal recorded by the pixels, it is possible to extract the g parameter. The advantage of this technique is that a g 
parameter can be estimated for each FD site, and it can cover lower scattering angles as shown in Fig. [6] (right) (the 
angular range where larger aerosols could be detected). It gives a similar estimation for the asymmetry parameter: 
(g) =0.54 ±0.16 (see Fig. [6] (left)). 



3.3 Wavelength dependence 

Up to now, the measurements to estimate the aerosol optical depth have been done at a fixed wavelength, that of 
the laser of the CLF. The empirical parameterization defined in Eq. ([3| gives a basic way to get the aerosol effect 
for any wavelength if the Angstrom exponent is known. At the Pierre Auger Observatory, two different facilities were 
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Fig. 6. Aerosol size measurements by the APF monitors and the very inclined CLF shots. Left panel: Asymmetry 
parameter distribution measured by the APF monitor (continuous line) and the very inclined CLF laser shots (dotted line). The 
APF monitor data were measured at Coihueco between June 2006 and July 2008 |5J. The values from the new method were 
obtained with laser shots at Los Leones recorded in 2008 [23]. Right panel: Relationship between the asymmetry parameter of 
the modified Henyey-Greenstein function and the mean aerosol size, i.e. the diameter in the case of spherical aerosols. Upper 
limits on the aerosol size are plotted for the two methods, following the Ramsauer approach |17j . 
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Fig. 7. Normalised distributions of the Angstrom exponent observed by the two monitors dedicated to this 
measurement at the Pierre Auger Observatory. Estimations of 7 from the HAM for data recorded between July 2006 
and February 2007 (continuous line) [13] ■ The measured 7 distribution for data collected by the FRAM from June 2006 until 
December 2008 (dotted line) [14]. 
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Fig. 8. Aerosol sampling at the Coihueco site (from |24J). Left panel: Andersen-Graseby 240 dichotomous sampler 
located at Coihueco, on the roof of the FD building. Right panel: Aerosol mass concentration for PMio measurements at the 
Coihueco sampling site (the circles). The horizontal lines give the mean values for the three aerosol samplings. 

installed to estimate this key parameter: the Horizontal Attenuation Monitor (HAM) and the F/ (Ph)otometric Robotic 
Atmospheric Monitor (FRAM). 

The HAM consists of a high intensity discharge lamp located at Coihueco, providing a large wavelength range 
source. To record the emitted light, a CCD camera is placed around 45 km away at Los Leones |13j . With this 
configuration, the total horizontal atmospheric attenuation across the Auger array can be obtained. Using a filter 
wheel, the camera records the aerosol attenuation between the two sites at five different wavelengths between 350 and 
550 nm. By fitting the signal with respect to the wavelength, it is possible to estimate the Angstrom exponent 7. From 
the data recorded between July 2006 and February 2007, the wavelength dependence displays a value 0.2 < 7 < 1.2 
(see Fig. m). The uncertainties are dominated by measurement fluctuations, and include a systematic effect due to 
subtraction of the estimated molecular attenuation between Los Leones and Coihueco. 

The main task of the FRAM telescope is the continuous monitoring of the wavelength dependence of the total 
column aerosol optical depth [T3]. The facility is located close to the Los Leones FD building. The method is the 
following: the telescope observes a set of chosen standard stars and, from these observations, it obtains extinction 
coefficients and their wavelength dependence. Five filters are used, for wavelengths between 360 and 547 nm. As 
preliminary results, the FRAM telescope estimates a value —0.5 < 7 < 0.3 (see Fig. [7]), lower than the results from 
the HAM instrument. 

The small value of the exponent 7 suggests that the Auger Observatory has a large component of coarse-mode 
aerosols, meaning aerosols larger than around 1 /um. 



3.4 Aerosol sampling measurements 

So far, all the techniques mentioned before provide information on aerosols as attenuators or scattering centres, but do 
not characterise the aerosols themselves. Such a characterization could help in the understanding of aerosol behaviour 
in the attenuation process. For instance, the atmospheric aerosols at the Pierre Auger Observatory are assumed to be 
desert-type particles [13] • Also, in the air shower reconstruction done in the Auger Offline software, the aerosols are 
pure scatterers: absorption is not taken into account, since it is assumed to be very small. But this assumption depends 
on the chemical composition of the aerosols. The results from aerosol sampling add new information about the aerosols 
at the site and can be compared with measurements obtained by the instruments using optical techniques such as the 
CLF or the lidars. They can give a clearer idea about the origin of the aerosols present at the Auger Observatory, their 
sources and trajectories (see Sect. [4]), and the connection between their composition and the meteorological variables in 
the region. An Andersen-Graseby 240 dichotomous sampler (see Fig. [8] (left)) provided with polycarbonate membrane 
filters was used to separate fine size particles PM2.5 (i.e. with an aerodynamic^ diameter d < 2.5 /im) and coarse 
particles PM10 (d < 10 /im) [24 . The sampling period was one day, beginning at 12:00 a.m. and ending at 12:00 a.m. 



1 Aerodynamic diameter is a physical property of a particle in a viscous fluid such as air. In general, particles have irregular 
shapes with actual geometric diameters that are difficult to measure. Aerodynamic diameter is an expression of a particle 
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Fig. 9. Micrograph using Scanning Electron Microscopy (from [12j). Left panel: Sampling of PM2.5, collected on 7 th 
July 2008. Right panel: Sampling of PM2.5-10, collected on 27 th October 2008. 
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Table 1. Element concentrations of PM2.5 and PM2.5-10 fractions, obtained by SEM/EDX (from [24J). Average 
particle concentrations, or atomic percentages, applied to two samples, one of each PM fraction. 



of the following day. The samples were collected from June to November 2008, at the Coihueco location, on the roof 
of the FD building. 

For gravimetric analysis and concentration measurements, a total of 34 filters containing fine particles (PM2.5) 
and 38 filters containing coarse particles (PM2.5-10) were considered. The concentrations, expressed as /zg/cm , were 
calculated as the ratio between the collected mass and the volume of air passed through the sampler during a full 
sampling period. Figure [8] (right) gives the evolution in time of the measured aerosol mass concentrations. An increase 
of the aerosol concentration is observed from winter to spring, which could be related to decreasing snowfall and 
increasing temperature. Although snowfalls are rare during winter, the low temperatures keep the snow from melting 
for long periods, reducing atmospheric particulates near the ground. 

Particle morphology and elemental composition were studied using a Scanning Electron Microscope (SEM) equipped 
with an Energy Dispersive X-ray system (EDX). Table [l] gives average compositions determined by SEM/EDX on 
at least 30 individual particles. It shows that Si, Al, Ca, Mg, Fe and K, the typical mineral soil elements, are the 
major components. This indicates that aerosols are mainly dust suspended in the atmosphere. PM2.5 fractions show 
more undetected particles due to the fact that beam focusing precision is insufficient to get a proper X-ray signal on 
particles smaller than one micron; thus many of them cannot be detected. Also, some particles in both fractions giving 
no detectable X-ray signal may be composed of light elements (Z < 11), presumably organic matter. SEM micrographs 
of two representative PM2.5 and PM 2 .5_io samples are shown in Fig. [9] 

4 Origin of atmospheric aerosols at Auger site using backward trajectory of air masses 

The Pierre Auger Collaboration has developed a large atmospheric monitoring program to have a better knowledge 
of aerosols present over its array. However, the aerosol population is highly mobile. Aerosol properties for a particular 
location are a function not only of local but also of emissions and meteorological conditions in the surrounding area. 
In our case, the Pierre Auger Observatory is affected by air masses potentially carrying aerosols from, for instance, 
Chile or the Pacific Ocean. Thus, any study having as a goal to explain aerosol properties over the Auger array needs 



aerodynamic behaviour as if it were a perfect sphere having a density of 1 g/cm 3 and diameter equal to the aerodynamic 
diameter. 
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Fig. 10. Monthly frequency all along the year of clear and hazy nights. Aerosol optical depths at 3.5 km above 
ground level measured between January 2004 and December 2010 at Los Morados. Clear nights are defined as 0.00 < r a < 0.01 
(continuous line) and hazy nights as r a > 0.10 (dotted line). 



to include work on the air mass origins. For a region outside of the Auger array to influence local measurements, the 
following three conditions must occur simultaneously: the region is a source of aerosols, the air mass coming from the 
region has to transport the aerosols in a sufficiently short time so that the aerosol content is still present on arrival, 
and the air mass arriving at the Auger array must be at an altitude that would affect the part of the atmosphere 
probed by the detector. Meteorological conditions also affect the aerosol transport to the observatory location. 

Different air models have been developed to study air mass relationships between two regions. Among them, the 
HYbrid Single- Particle Lagrangian Integrated Trajectory model, or HYSPL1T 25 , 26J , is a commonly used air- modelling 
program in atmospheric sciences that can calculate air mass displacements from one region to another. The HYSPLIT 
model developed by the Air Resources Laboratory, NOAA^J is a complete system designed to support a wide range 
of simulations related to regional or long-range transport and dispersion of airborne particles. HYSPLIT computes 
simple trajectories to complex dispersion and deposition simulations using either puff or particle approaches with a 
Lagrangian framework. HYSPLIT can be used to get backward/forward trajectories: by moving backward/forward 
in time, the resulting backward/forward trajectory indicates air mass arriving at a receptor for a particular time, 
identifying the regions linked to it. 

From values of aerosol optical depth given in Fig. [4] at Los Morados, the data sample is divided into three popu- 
lations: the clear nights with the lowest aerosol concentrations (r a < 0.01), the hazy nights with the highest aerosol 
concentrations (r a > 0.10), the average nights being in the third category. Figure 10 shows the relative frequency 



month-by-month for clear nights and hazy nights. Clear nights are more common during the Austral winter than the 
rest of the year. Using the HYSPLIT tool, each clear or hazy night is assigned to a backward trajectory computed 
over 48 hours. Meteorological quantities used in calculations come from the GDAS model (Global Data Assimilation 
System), its validity having already been checked by the Pierre Auger Collaboration [2"T] . 



2 NOAA, National Oceanic and Atmospheric Administration, U.S.A. 
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Distribution of trajectories: 0.00 < i a (3.5 km) < 0.01 Distribution of trajectories: 0.10 < x a (3.5 km) < 0.30 
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Air mass path directions: 0.00 < x a (3.5 km) < 0.01 Air mass path directions: 0.1 < T a (3.5 km) < 0.30 




Fig. 11. Distribution of backward trajectories and direction of air masses for clear nights and hazy nights 
(from [28j ) . Paths estimated with HYSPLIT for the years 2007/2008/2009 and aerosol optical depth data coming from the 
CLF measurements. Top panel: Distribution of 48-h backward trajectories from the Malargiie location for a start altitude fixed 
at 500 m AGL. Bottom panel: Direction of trajectories at three different start altitudes AGL: 500 m [solid line in blue), 1000 m 
[dashed line in red) and 3000 m (dotted line in magenta). Each distribution is normalised to one. Colours in online version. 



The distributions of the backward paths for clear nights and hazy nights are shown in Fig. [TTJ (top). The two 
plots are quite different: the air masses come mainly from the Pacific Ocean during the clear nights and travel 
principally through continental areas during the previous 48 hours for hazy nights. Following the conclusion of the 
chemical composition aerosol analysis presented in Sect. |3.4[ the aerosols originate mainly from the soil. Thus 48-h 
backward trajectories traveling mainly over the ocean can be characterised as air masses with a low concentration 
of soil aerosols. The aerosol contamination occurs only during the last hours for "oceanic" trajectories, whereas this 
is possible during the full period for the overland trajectories. It seems that backward trajectories are a key point, 
explaining the differences in aerosol concentration at the Pierre Auger Observatory. In Fig. [TTJ (bottom) , directions of 
air masses are given for clear and hazy nights, for three different start altitudes at the location of the Pierre Auger 
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Fig. 12. Example of a forward trajectory from the Malargiie location using HYSPLIT. The initial height is fixed 
at 500 m above ground level and the time scale is 100 hours (map taken from Google Earth). Kerguelen Island is represented 
with the circle. 



Observatory: 500 m AGL (blue diagram) , 1000 m AGL (red diagram) and 3000 m AGL (magenta diagram) . The same 
conclusions can be drawn for all three initial altitudes. Even if a slight shift to the West is seen for the highest altitude 
for hazy nights, it seems important to remind oneself that atmospheric aerosols are usually located in the low part of 
the atmosphere, typically in the first 2 km. 



5 The Pierre Auger Observatory as an aerosol observatory for Southern Hemisphere 

We present in this section an interdisciplinary proposal for using the Pierre Auger Observatory as a scientific facility 
that will considerably enhance our knowledge of aerosol cycles over Southern Hemisphere. It has been shown in this 
text that the Pierre Auger Collaboration has developed a large atmospheric monitoring program with, currently, 
almost seven years of accumulated data [35] ■ Beside the benefit for astroparticle measurements, a better knowledge of 
aerosol properties will enhance our capability of climate change forecasting because of the impact of dust deposition on 
oceanic biogeochemistry. We would obtain a better accuracy on atmospheric fluorescence measurements plus valuable 
new data concerning continental dust emission, transport and deposition. The next sections explain the scientific issues 
of the proposed work. 



5.1 Impact of aerosol deposition on the oceans 

Atmospheric mineral aerosol, so-called "dust" , is one of the major vectors feeding open ocean surface waters with 
trace metals [30 . Even at extremely low concentrations, trace metals are micro- nutrients, necessary for the growth of 
phytoplankton [31] . In this way, trace metals are linked to climate change since they affect the capability of the marine 
biomass to trap CO2 (see for example [32] )• The Austral region ranging from about 40° and 65° South is a major 
oceanic CO2 sinks with exported carbon to deep sea water |33j . This region is also very remote from continents and thus 
atmospheric dust exhibits very low concentrations [34,35,36,37]. This oceanic area is a HNLC region (High-Nutrients 
Low-Chlorophyll) and dust deposition could be a severe limiting factor for the primary production of nutrients [3H1 
[39] . Deposition measurements of continental dust have never been performed in this large area, except in the 1970's 
associated with radon [40 . Some field measurements are available from point measurements at Crozet Island [4"T] . 
Others use the deposition flux calculated from aerosol concentration measurements using rough assumptions based 
on ad-hoc deposition velocities or scavenging ratios |37j . Recently, it has been shown by performing direct deposition 
measurements on Kerguelen Island that this last approach leads to atmospheric fluxes 50 to 100 times lower than 
those directly measured |42j . There is an on-going debate on the importance of the atmospheric deposition flux 
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over the Southern Ocean. Early modelling studies hypothesized a large flux, whereas in- situ measurements have led 
to a downward revision these early estimates. More recently, direct flux measurements again have pushed up the 
overall dust deposition flux to the Southern Ocean. Cassar et al. [38] stated from modelling studies that the primary 
productivity is proportional to the aeolian dissolved iron deposition originated from dust. Iron bioavailability, and 
those of other nutrients, is conditioned by the ability of mineral dust to be dissolved in seawater. In the current 
estimation of atmospheric dissolved iron fluxes, the solubility of atmospheric iron in seawater remains a major source 
of uncertainty [43] . Fluxes are generally calculated with an average solubility value obtained from available data in 
literature [371138] . However the large majority of data are related only to African and Asian dust. Given the key role 
of the Austral Ocean on global climate (as a main oceanic CO2 sinks), it is important to do studies that aimed to 
characterize iron and other trace metal solubility from Southern Hemisphere dust. 

Since the Southern Hemisphere is mainly constituted of oceans, the only possible dust sources are South America 
including East Patagonia and East Andes Foothills, South Africa and Australia. South America is suspected to be the 
major dust source for the oceanic region ranging between 40° S and 60° S, with a major contribution of Patagonia but 
also with important sources located more North in the arid or semi-arid regions located East of the Andes mountains. 



5.2 Regional emission of mineral aerosol (dust) 

The Pierre Auger Observatory may be very useful for environmental science in general and particularly those concerning 
aerosols. The observatory is located in an area subjected to local emissions of soil dust but also to long-range transport 
of aerosols having marine, terrigeneous or anthropogenic origins. Such transported aerosols are often located at high 
altitudes and can be detected by laser beam scattering, in-situ aeroplane, kite or balloon sampling. They reach ground 
level during deposition processes. On the other hand, locally emitted aerosols are more present at ground level and can 



be well characterized by performing ground-based sampling on filters (see Sect. 3.4). This should allow investigations of 
the chemical composition and particularly of the content in iron bearing particles that can absorb blue light [44] . These 
aerosols can reach higher altitudes and be forward transported far away in the general sub-Antarctic atmospheric 



circulation (see Fig. 12). Because of the sporadic character of aerosol emission pulses, continuous sampling should 
be carried out. Transportation models, including air mass trajectories, can provide a framework to understand the 
behaviour of the aerosol loading the atmosphere over the observatory, and assessing impacts on the Southern Ocean. 



5.3 Aerosol characterization and proposed action at the Pierre Auger Observatory 

In the previous sections, we have pointed out the influence of the aerosol and deposition chemistry on their potential 
impacts. Routine measurements are already being done using lidar systems, and particle number size distribution 
determined using optical particle counting. For mass balances, the particle mass size distribution will be deduced 
using model that makes assumptions about density and particle shape |3S]. Some preliminary sampling campaigns 
were done, using filtration or impaction. For filtration, air is pumped through a filter, often a membrane with a 
0.2 — 0.4 /im porosity. For impaction, a narrow air jet is directed against a plate, the overall geometry of the system 
will separate by impaction particles following their aerodynamical diameter. The larger particles impacting on the 
plate stay attached, whilst the smaller ones do not. If several of such impactor plates are cascaded, we obtain a so 
called cascade impactor which can deliver up to 13 size segregated levels of particles with a decreasing aerodynamical 
diameter. With a simple computation method, a continuous size distribution in mass can be rebuilt [46] . Filters and 
plates must be brought to laboratory where various chemical analyses are performed. Because of the high variability of 
atmospheric conditions and compositions, we propose to set up a continuous time series of aerosol sampling, allowing 
a quantitative annual budget. This continuous aerosol survey will run at Coihueco station, on a weekly basis. Various 
measurements will be done on filters: elemental analyses and more complex chemical analyses, including solubility. The 
only limitation is the preservation of the sample, which is not an issue except of contamination. In-situ experiments 
will use only small pumping systems and optical counters and do not need more than some hundreds of watts of 
electrical power. Most of time, work will be done at the laboratory. 

For elemental analyses of collected samples, thin layer X-Ray fluorescence spectrometry is operated in laboratory 
without any transformation of the sample 47,45], and at high rate and low cost. PIXE (Proton Induced X-Ray 
Emission) can also be used [4"8"ll4T)] . To access ultra-trace amounts and isotopic determinations, samples have to be 
digested in a mixture of various acids, to put all the sample in a soluble form. Elemental analyses are performed by 
Inductively Coupled Plasma Atomic Emission Spectrometry (ICP-AES) [50] and Inductively Coupled Plasma Mass 
Spectrometry (ICP-MS) [EE]. 

Elemental analyses will be performed on all the collected samples for major and trace elements using XRF. These 
first analyses will be used to select interesting samples. On these samples, ultra-trace and solubility measurements, 
which require more work, will be done using ICP-AES and ICP-MS. For solubility, water at various pH and various 
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composition, including sea water, is passed through the sample and dissolves parts of the aerosols in a kinetic experi- 
ment |52U53j . Variations of dissolution rate along time but also acidity of the water give valuable information on the 
bioavailability of the collected aerosols but also chemical bonding in the solid aerosol. To enhance the interpretation of 
such chemical data, particles will be directly observed by Transmission Electronic Microscopy (TEM) and Scattered 
Electronic Microscopy (SEM). Microscopy is devoted to single particle analyses, and observation of numerous particles 
can give very detailed information on size and shape distribution, chemical composition and mineralogy of the collected 
aerosols |54j . Such information, together with bulk elemental analyses and solubility measurements are valuable for 
inclusion in a global model of dust solubility and dust impact on biogeochemical cycles |55j . Such knowledge will be 
finally used to enhance lidar and CLF measurements interpretations and to improve our knowledge of the influence of 
aerosol chemical composition on the astroparticles atmospheric signals. Indeed, for instance, it would be interesting 
to check if assumptions on aerosol absorption done up to now by the Pierre Auger Collaboration are validated. 



Conclusion 

The Pierre Auger Collaboration has accumulated a large database of aerosol measurements. This effort has signif- 
icantly reduced systematic uncertainties in the cosmic ray air shower reconstruction. Cross-checks between optical 
measurements and sampling, including their size distribution or their chemical composition, are now being undertaken 
at the observatory to improve our knowledge on aerosols. 

Collaborative works with scientists from the environmental sciences are currently discussed in the Pierre Auger 
Collaboration. Indeed, the observatory is located in a region where only few atmospheric survey stations are installed. 
With additional facilities as enumerated in Sect.[5j the Pierre Auger Observatory could contribute to a better knowledge 
of the Southern Hemisphere atmosphere and its chemical composition. The observatory could be candidate to become 
member of international networks as the Global Atmosphere Watch (GAW) program [S3] or the Aerosol Robotic 
Network (AERONET) [57] and more generally to be a reference station. 

The proposed collaborative project would considerably enhance knowledge of the optical and chemical properties of 
aerosols over the Pierre Auger Observatory but also the global impact of airborne particles over Southern Hemisphere. 
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